Pure spin decoherence in quantum dots. 
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We uncover two microscopic physical settings with significant pure spin decoherence. First, for 
quantum dots (QD) electrostatically confined in two-dimensional hole gas, decoherence comes from 
qubit spin-orbit (SO) coupling to phonons, whose decay due to free charge carriers in contacts 
is described by Ohmic weighted phonon spectral function. We derive significant SO interactions 
affecting holes with origin and symmetry distinct from that of conventionally considered Dresselhaus 
and Rashba terms. In the second setting of electron or hole QDs coupled to a linear chain of atoms, 
decoherence is due to spin-dependent coupling to phonons, whose decay due to scattering off the 
free ends of a chain is described by the weighted phonon spectral function inverse proportional to 
frequency. The decoherence rate in both settings is linear in temperature. 

PACS numbers: 73.21.La, 71.70.Ej, 03.65.Yz, 03.67.-a 



When spin population is conserved, quantum coher- 
ent properties are defined by pure spin decoherence pQ. 
In optics, there is a related phenomenon of zero phonon 
spectral lines [2H3]. For the field of quantum computing, 
understanding mechanisms of the pure spin decoherence 
is vital. In quantum dot (QD) qubits, the task is to un- 
cover qubit-phonon interactions and processes in phonon 
system that lead to the effect. Current understanding 
of QDs does not include pure spin decoherence due to 
spin-orbit (SO) interactions and phonons as an effective 
mechanism, and finding to the contrary is important. 

In this letter I show that pure spin decoherence is cru- 
cial for QDs. In particular, I consider QD qubits with 
charge carrier holes, in which spin relaxation is sup- 
pressed, and nuclear spin effects are weak compared to 
electron QDs. I derive spin-phonon interactions in hole 
QDs giving sizable pure spin decoherence. I show that 
steps can be taken to suppress it. I find that pure spin de- 
coherence in QDs can be related to realistic microscopic 
mechanisms of phonon decay. It arises for a hole localized 
in a QD electrostatically confined in the two-dimensional 
(2D) gas, when the hole interacts with phonons, which 
in turn decay due to interaction with free holes present 
in the system, Fig[l] This is a microscopic realization 
of models with Ohmic phonon weighted spectral den- 
sity. The other setting for pure decoherence, Fig[2] is 
electronic or hole QD coupled to ID atomic chain, due 
to qubit-phonon coupling and phonon scatttering of the 
free ends of the chain, which results in a weighted spec- 
tral density inverse proportional to frequency. 

I derive the effective Hamiltonian for QDs confined to 
2D hole systems, which also describes rings, wires and 
point contacts. The problem was extensively addressed 
in recent years [5Tll6|. but an important effect was missed. 

In the ground state of the confined hole systems, like 
heterostructures in iii-v or Si/Ge structures grown along 
001 crystallographic direction, the effective mass rrih in 
the growth direction z is a heavy hole mass. This state is 



FIG. 1: QD in a 2D hole gas. Pure spin decoherence is due to 
interaction of a localized hole spin with phonons that in turn 
decay due to interaction with delocalized holes in contacts. 
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FIG. 2: a. QD coupled to a chain of atoms. Pure decoherence 
is due to spin coupling to phonons that scatter of the free ends 
(0,1) of a chain, b. Inductance-capacitance circuit equivalent 
to an ideal chain, c. Circuit for a chain with missing spring. 

attractive for qubits [5], because it is mostly character- 
ized by the angular momentum projection L z = ±3/2. 
Since the matrix elements of spin-flipping interactions 
due to J x and J y between states +3/2 and -3/2 vanish, 
decoherence of hole QDs is reduced [TTJ [TH] compared 
to electronic QDs. Understanding hole qubits requires 
knowledge of the effective Hamiltonian projected on the 
hole ground state doublet that must account for mixing 
of ±3/2 states to the states with light effective mass m; 
in z-direction ( L z = ±1/2). 

To date, SO interactions for 2D and QD holes have 
been treated similar to those for electrons, and limited 
to modified Dresselhaus and Rashba terms, which do not 
give pure spin decoherence. Here we show that SO inter- 
actions distinct from the Dresselhaus and Rashba terms 
affect 2D and QD holes. They arise even if bulk crystal is 
ccntrosymmetric (no Dresselhaus term |19j). and quan- 
tum wells (QW) are symmetric (no Rashba term j2"0"]). 
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Their origin is the hole Luttinger Hamiltonian 



H = 



2m 
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where Ji are the 4x4matrices of the i— component of 
J = 3/2 operator, [AB] = (AB+BA)/2, 7i are Luttintger 
constants, mo is the free electron mass. Infinite rectangu- 
lar QW, in which Hamiltonian Eq.Q is supplemented by 
the boundary conditions giving vanishing wavefunctions 
at z = ±rf/2, allows analytical solution 22-24]. Wave- 
functions ^s(z), ^a(z) are symmetric and antisymmetric 
with respect to reflection about the plane z = 0, and in 
the basis of Bloch functions (m 3 / 2 , u 1 / 2 , w -1 / 2 , u~ 3 / 2 ) de- 
scribing bulk holes with J = 3/2, are given by ^± (z, r) = 



a± z exp (ik • r), where r is the in-plane coordinate, and 
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k is the in-planc momentum, 4>k is the angle between 
k and x— direction. Here S z = sing^z — (sh/si) sin qi z, 
C z — cosq^z — (ch/ci) cos qiz, are antisymmetric and 
symmetric standing waves with respect to z — 0, = 
sinq/jd/2, sj = smg;<i/2, Ch — cosqhd/2 c/ = cosqid/2, 
and d is the QW width. For energies of in-plane mo- 
tion ek < e g = h 2 q 2 /2rrih, the wavevectors qh = ir/d 
and qi = y/vqh, where v = mi J rah,. Thus, QW holes 
in contrast to QW electrons are described by two stand- 
ing waves, illustrating mutual transformation of heavy 
and light holes upon reflection from the barrier. In the 
spherical approximation 72 = 73 and at small in-plane 

energies, the normalization factor A = I ~ + 4^/ s 1 , 
/. = (l+^)f + ft^§, and amplitudes ^ =A #, 



A 2 = A ^gf and A 3 = 3A k 3 /{8vq 3 h ). 

Taking the projection [35] of V(r, z, t) that includes po- 
tential confining holes to QDs and potential due to long- 
range hole-phonon interactions 26, .27] onto the ground 
state given by Eqs.Q, in the lowest order in k we derive 

n = ts£ + "oVb(r, t) + aiVi(r, t) - ia 2 k • V r Vi(r, t) 
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a 2 a z [k x V T Vi(r,t)] z , 
where m* is the effective mass given by 
m 72(71 + 72) - H , 3 73 2 (7i 2 - 47I) 1 / 2 
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where toq is the free electron mass, f{6) = ^ , 6 
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jection potentials Vb(r, t) = | J^f, 2 V((r, z,t)C 2 (z)dz, 



Vi(r,t) = lf* / d 2 /2 V((r ) z,t)S 2 (z)dz. The SO terms in 
([3]) give no spin-flips, and have U(l) symmetry, partially 
breaking spin rotation symmetry [28j . Spin flip terms 
from projection of V(r, z, t) onto ^± (z, r) contain at least 
the parameter (kd) 3 and are smaller than SO term in |3j 
at small in-plane energies. We note that for 2D holes, the 
Rashba term is oc fc 3 , and the Dresselhaus terms include 
both cubic and linear in k contributions [29]. Thus, the 
SO term in ^ missed in all previous work on hole QWs 
and QDs has profound consequences: in Si-Ge systems, 
it dominates; in iii-v systems, it changes the symmetry of 
SO interactions. This term is of sizable strength, defined 
by kd and averaging over z-direction, and is not related 
to admixture of distant bands like electron SO terms. 

We now discuss pure decoherence in small QDs, with 
sizes less than hole mean free path. We take parabolic 
projection potential aoVo+aiVx — m* (uj 2 x 2 +uj 2 y 2 ) /2 to 
model in-plane spatial confinement of holes. The linear 
in k and V% term gives negligible mixture of orbital levels 
and is omitted. The relevant SO interaction is the last 
term in |3]) due to hole-phonon interaction V(z, r, t). We 
will refer to it as spin-phonon interaction 



E 
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where Vqj is the Fourier component of the amplitude of 
spin-phonon interaction, j denote the phonon mode po- 
larization, and c q j , c q j are phonon creation/annihilation 
operators, which define lattice displacement operator 
Rnp = S q ,j u « P (q, j)(c qj - + C-qj), where n labels the 
crystalline unit cell, p are the atoms in the cell, and 



un P (q,i) = 



2 P Vuj, 



1/2 



A .piq-x(np) 

u qj e > 



(6) 



is the amplitude of the phonon in a mode q, j at a site 
(np), V is the crystal volume, d q j is the phonon polar- 
ization vector, and x(np) is the lattice vector. 

In calculating decoherence we assume the presence 
of uniform magnetic field H||z with Larmor frequency 
u) z . The projection of the spin operator a z is conserved. 
Hence, only pure decoherence can affect the qubit, once 
it is in a superposition of up and down state, with the 
pure decoherence constant r_|_ defining the motion of the 
transverse average spin: 



S x — (w z x S) x — T±S X . 



(7) 



For the spin-phonon interaction ([5]), Tj_ = T xx — T yy . 

Our task is to find a microscopic mechanism(s), which 
result in non-zero T±, and to calculate this constant. It is 
important to realize that spin-phonon interaction |5]) will 
not lead to pure decoherence on its own. The effective 
magnetic field that acts on transverse spin due to this in- 
teraction is time and coordinate-dependent. If phonons 
are coherent and do not decay, they would result in co- 
herent modulation of the ground and excited QD Zeeman 
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energy levels, with time-dependent phase. Indeed, at a 
given point in space, spin interaction with phonons is pe- 
riodic in time. One frequent approach to decoherence is 
to introduce a phenomenological decay/correlation time. 
As we shall see, the overwhelming majority of phonon de- 
cay processes do not give T±. In particular, although it 
was suggested [3U] that spin-phonon interactions give de- 
coherence because bulk phonons enter and leave the QD, 
such process cannot contribute on its own. The fact that 
spin-phonon interaction occurs only in a QD is accounted 
by the qubit wavefunctions defining Vq.j , and additional 
factors of phonon decay need to be considered. 

To find out limitations on processes leading to pure 
decoherence, we calculate T±, which is defined (see, e.g., 
PQ) by the Fourier component of the correlator of fluctu- 
ating spin-orbit fields at zero frequency 



T x (uj = 0) 



(n,(0)n,(r))dr, 



(8) 



where angular brackets denote statistical average. The 
correlator in Eq.Q is defined by the Fourier im- 
age of the phonon correlation function C(ui) — 
§E„ |^ q | 2 ( Cq (r)c_ q (0))+h.c., given bygH E2] 



C(«) 



(2N(io) + 1)1%) 

( W 2- Wq )2/ w 2 +r p( w )- 
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Here N(oj) is the number of phonons with energy tuo, 
and n(cj) is the weighted phonon spectral function (in 
terms of Keldysh formalism, the imaginary part of re- 
tarded phonon self-energy). Real part of phonon self- 
energy, renormalizing phonon energy, does not result in 
decohehence and is neglected. We aim at finding H(uj) on 
rigorous microscopic grounds, without introducing any 
phenomenological parameters of phonon decay. 

The physics of dependence of T± on Il(w = 0) is 
that no energy exchange between Zeemann levels and 
phonons occur. SO hole-phonon interactions are only 
shaking up Zeeman levels. Physics is reminiscent of 
that of the zero-phonon line width in exciton or molec- 
ular absorption/emission [5H1]. However, spin dephas- 
ing and exciton linewidth problems differ. In the lat- 
ter, carrier-phonon interactions give linear and quadratic 
in the phonon coordinates terms in two-level systems, 
and quadratic terms have been considered the princi- 
pal cause of zero-phonon linewidth. In the former, the 
SO carrier-phonon term linear in phonon coordinates de- 
pends on er 2 , its square is spin- independent, and only the 
linear term can contribute to Tj_(uj = 0). Some works 
[31 0] conclude that interactions linear in Vq do not con- 
tribute to linewidth, or contribute only in the case of 
model localized phonons, rather than for realistic phys- 
ical phonon/ vibration modes. We now show two cases 
with real phonons contributing to relevant T±(lu = 0). 

In case I, n 1 ^) = w/(2w q r ( J, (i ). Here l/r^ is the 
phonon decay rate at u> — u; q for mechanisms giving 



1/t<J, oc io. This is the case of Ohmic phonon spectral 
function, discussed in the theory of dissipative quan- 
tum tuneling by Leggett et all [33] . Here we calculate 
= 0) due to the process, in which phonons that 
interact with QD spin decay as a result of the phonon 
scattering off delocalized holes neccesarily present in the 
system with electrostatically confined QDs. The case II, 
in which C'(lo = 0) contributes to T±, was never treated, 
and occurs when H ii (uj) — w q /(2wr q I ) Here 1/t^ 1 is 
phonon decay rate for mechanisms giving oc ljuo. 

We find the mechanism due to phonons interacting with 
qubit coupled to ID linear chain, in which phonons decay 
due to scattering of the free ends, giving the case II. 

By definition, Tl{uS) is an odd function of lo. It is easy to 
see that cases I and II give the only possible dependencies 
of n(cj) that contribute to non-zero meaningful C(u> — 
0) = lim^^oC (uj) . The resulting T±(uj = 0) is 
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case I, 
case II, 



(10) 
(11) 



T is the temperature, kg is the Boltzmann constant. Re- 
markably, T± oc T. This contrasts with known mecha- 
nisms of spin relaxation and decoherence in QDs, but is 
not surprising for pure dephasing, which, at uj = 0, is 
similar to classical mean field fluctuations [34] . We note 
in passing that spin-independent QD electron-phonon in- 
teraction with consequent phonon decay due to electrons 
in contacts gives case I pure decoherence that can explain 
T^ 1 dependence of dephasing time in experiments [35]. 

We now turn to microscopic calculation of T± . 
Case I: In GaAs systems, the piezoelectric interaction of 
holes with acoustic phonons that defines the last term in 
Eq.((3j and spin-phonon interaction Eq. ^ is given by 



V(v,t) 



E 

qj 



h 



2pVu, 



-h.c. 



1,3 



(12) 

where = /9|e<jy| is a piezoelectric constant, Ej/y is 
an antisymmetric tensor, ftcaAs = 1.2 X 10 7 eV/cm, = 
qi/q. Denoting A; = y // h/m*u>i the sizes of the wavefunc- 
tions <f> in ID parabolic potentials, we find the QD orbital 
ground state wavefunction "J ~ ^oo+^iij where ^ nm = 
<t>n{x)4>m{y), and a = iuj c {\ 1 /\ 2 - A 2 /Ai)/4(cJi +u 2 ). In 
weak magnetic field, the lowest QD states are due to Zee- 
man doublet of the ground orbital state. 

The hole-phonon interaction V% in Eq.|3~| defined by 
( 12 ) gives 2 oc 



(y*\e^{q x ky-k x q y )\y) = a{q 2 x X x /\ y - 
qyX y /X x ). It also leads to factor A(q) = |d qj | 2 en- 
tering | Vq | 2 in Eq. (|9|). In a simple model with isotropic 



speed of sound s, A(q) = {q 2 q. 



j 2 



)M 



Calculating n(w) due to interactions of phonons with 
holes in contacts, we neglect renormalization of the hole- 
phonon vertex, like in consideration |31) of electron- 
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phonon interactions in 3D metals. In the 2D case 

n» = J G(p ) e)G(p-q,e-u } )d 2 pde (13) 

where G(p, e), ep and pp are the Green function, Fermi 
energy and momentum for holes in contacts. At small u) 



from the solution of Dyson equation [JO] we find 



/3 2 rn*A(q) p F Tiw 
Anhdpujq q ep 



Q(2p F -q), (14) 



where <d(x) = 1 at x > and ®(x) = at x < 0. 

Calculating the hole ground state diagonal matrix el- 



ement of (12), evaluating Eq.(14), and taking ui — > 
in Eq.([9|, for to* = 0.25too, hole density in contacts 
p = 6 x 10 11 cm" 2 , s = 4 x 10 5 cm/s, d = 100A, 
\ x = 400A, X y = 250A, i? = 0.5T and T = 100m.fr, 
from Eq.rfToJ we get r ± (w = 0) = 10 4 Hz. Thus, pure 



decoherence is significant in GaAs QD. This rate ex- 
ceeds 10 times the hole decoherence rate experimentally 
observed in InAs QD To fight the pure decoher- 

ence, one makes QD shape close to circular. We note 
that in SiGe/Ge/SiGe QD structure, our estimates show 
r_i_(w = 0) = 10 3 Hz for similar QD sizes. With deco- 
herence due to nuclei in such structures weaker than in 
GaAs, pure dephasing proposed here becomes important. 
We address hole qubits in Si/Ge structures elsewhere [18]. 

Phonon decay due to anharmonic processes, scattering 
by impurities, surfaces, or QD or QW boundaries gives 
n(a>) oc uj s with s > 3 and does not lead to case I pure 
dephasing. Therefore observation of pure dephasing will 
make it possible to separate phonon decay due to charge 
carriers in contacts from all other decay mechanisms. 

Case II: This setting can be realized in QD coupled 
to ID linear atomic chain, where phonon decay occurs 
due to scattering off its free ends. Spin-phonon term 
([5]) can arise in both electron and hole QDs due to g- 
factor fluctuations [37]. For case II in electron QDs, 
we consider g-factor modulation in a QD of wavefunc- 
tion extent A by phonons of a chain of N atoms of 
mass M with lattice constant a linked by nearest neigh- 
bor springs of spring constant rj. Following |38j . we 

get interaction {5} V q = i{h/2MNuj q ) 1 / 2 Aqe- q2d2 / i hujz 1 
where ui 2 = (4r)/M)sin 2 (qa/2), and A is a constant, 
which can be obtained from spin resonance experi- 
ments [35]. We consider ends as defects in crystal with 
cyclic boundary conditions. In this model, 11* (u>) = 
J £q' Kuq'I^V + i0)|u q )| 2 <5(?iw - faj q ,). Here u q is 
given by Eq.6 in ID case, scattering matrix T(u>) with 
dimensionality of a spring constant is defined by Dyson 
equation T = SL + SLGT, G is the ideal chain phonon 
Green function, SL is perturbation in spring constants 
due to chain ends at atoms n = 1 and n = 0, with 
zero 7 of spring joining them. The nonzero SL nn i are 
SLqq = SL11 = —SLiq = —SLqi = T]. Expanding T(lo) in 
u q as T nn i(u) = 2M 2 /hJ2 quq2 u n (qi)t{qi,q2,uj)u n >(q2), 



t(q 1 ,q 2 ,uj 2 + £0) 



4?ye 



a (gi+Q2) 



(15) 



Eq.(15) gives t oc 1/ui that is key to singular weighted 



phonon spectral function and r¥: 



n^(w) = 2r]\smqa/2\/Moj (16) 
^ujzIujl) 2 ■ (kT/huL) ■ (A 2 h/MXa), (17) 



where u>\ = Arj/M . Physics of t oc l/u> resulting in such 
Uq 1 is clear from the electrical circuit equivalent to an 
atomic chain with missing spring, Fig 2c. At capacitance 
C = 00 (due to 7 = between ends and 1), the circuit 
is shorted, and u> — > responce to finite current J is 
zero voltage V (the responce to finite V is J oc 1/w). 
In a chain with free ends, finite force gives displacement 
oc 1/oj; when phonons scatter, t oc 1/ui. Taking A = 50 A 
for a carbon QD coupled to a chain of carbon atoms 
with M = 2 x 10~ 23 g, a = 2.5A, uj l = 1.2 x 10 11 Hz, 
A = 0.1, H = 0.5T, T = 100 mK, we get = 10 7 Hz. 
Such decoherence must be avoided. This shows, e.g., the 
limits on reducing bath dimensions in engineering qubits. 
Once phonons are 3D, case II is not realized. We note 
here that decay rates and are not additive: once 
case II decoherence is present, solely defines T±. 

In conclusion, we demonstrated that pure dephasing 
can considerably affect qubits in important settings. We 
presented two physical mechanisms based entirely on re- 
alistic microscopic quantum processes, without introduc- 
tion of any kind of phenomenological correlation time. 
We also demonstrated that important spin-orbit inter- 
actions distinct in origin and symmetry from known 
Rashba and Dresselhaus terms describe physics of low- 
dimensional holes. Support by NSF grant ECCS-0901754 
is gratefully acknowledged. 
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